symbol = {'1','2','3','4','5','6','7','8','9','*','0','#'};
lfg = [697 770 852 941]; % Low frequency group
hfg = [1209 1336 1477];  % High frequency group
f  = [];
freqs = [lfg hfg];

for c=1:4,
    for r=1:3,
        f = [ f [lfg(c);hfg(r)] ];
    end
end

FS=8000;
X=0:1/FS:0.1;
sigma2=0.1;
sigma=sqrt(sigma2);

RMSsin=1/sqrt(2);
SNR=20*log10(RMSsin/sigma2);

inputNoNoise = tone(2);
%input = awgn(input, (1/sqrt(2))/0.2);
input = awgn(inputNoNoise, -17);

[b1,a1] = ellipfilter(freqs(1));
[b2,a2] = ellipfilter(freqs(2));
[b3,a3] = ellipfilter(freqs(3));
[b4,a4] = ellipfilter(freqs(4));
[b5,a5] = ellipfilter(freqs(5));
[b6,a6] = ellipfilter(freqs(6));
[b7,a7] = ellipfilter(freqs(7));

out1 = filter(b1, a1, input);
out2 = filter(b2, a2, input);
out3 = filter(b3, a3, input);
out4 = filter(b4, a4, input);
out5 = filter(b5, a5, input);
out6 = filter(b6, a6, input);
out7 = filter(b7, a7, input);

out = [sum(out1.^2), sum(out2.^2), sum(out3.^2), sum(out4.^2), sum(out5.^2), sum(out6.^2), sum(out7.^2)]

